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On base of Hamiltonian formalism, we show that Hopf bifurcation arrives, 
in the course of the system evolution, at creation of revolving region of the 
C^ . phase plane being bounded by limit cycle. A revolving phase plane with a 

^ I set of limit cycles is presented in analogy with revolving vessel containing 

— Jj ' superfluid He^. Within such a representation, fast varying angle is shown to 

P^ ■ be reduced to phase of complex order parameter whose module squared plays 

O I a role of action. Respectively, vector potential of conjugate field is reduced 

to relative velocity of movement of the limit cycle interior with respect to its 
1 ■ exterior. 
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^ ■ 1 Introduction 

-4— > 

c^ I It is not gross exaggeration to assert that conception of the phase transition is one 
of fundamental ideas of contemporary physics. Related picture is based on the 
i-^ ■ Landau scheme, according to which a thermodynamic system, driven by slow and 

C , monotonic variations of state parameters of heat bath, rebuilds its macroscopic state 
if thermodynamic potential gets one or more additional minima in a state space P] . 
From mathematical point of view, such a phase transition represents the simplest 
bifurcation that results in doubling thermodynamic steady states. 

As is known, thermodynamic phase transition is a special case of self-organization 
process in the course of which three principle parameters, being order parameter, its 
conjugate field and control parameter, very in self-consistent manner |2]. Roughly 
speaking, a generalization of thermodynamic picture due to passage to synergetic one 
is stipulated by extension of set of state parameters from single parameter to three 
ones, being above pointed out. It might hope that such a generalization allows one to 
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describe not only the simplest Landau-like bifurcation, but much more complicated 
Hopf one, when a limit cycle is created to be continuous manifold instead of discrete 
one |Hj. 

Our considerations of this problem have shown 0] that using the whole set of 
universal deformations within standard synergetic scheme does not arrive at stable 
limit cycle, whereas running out off the standard scheme of self-organization has 
allowed us to obtain the limit cycle shown in Fig.l [S]. In this connection, the 
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Figure 1: The limit cycle related to non-standard self-organization equations x = 

x[y - (l + rx) -{A- l)/(l + px)], y = A-y{l + x) with A = 14, r = 5, p = 2 0. 

question arises: what is the physical reason that self-organization scheme may not 
represent the Hopf bifurcation? 

This paper is devoted to the answering above question. It is appeared, main 
reason is that a description of a limit cycle demands using both potential and force 
of a field conjugated to an order parameter, whereas standard synergetic scheme uses 
an force of this field only. Following [Hj, we show in Section 2 that fast revolving of 
the state point in the phase plane induces a gauge field, whose potential is reduced 
to relative velocity of movement of interior domain of the limit cycle with respect 
to the exterior. Such a picture allows us to study, in Section 3, a revolving phase 
plane with a set of limit cycles, using an analogy with revolving vessel containing 
superfluid He^ [2j. Section 4 concludes our consideration. 



2 Hopf bifurcation within canonical representation 

We consider Hamiltonian dynamics determined by the equations of motion 
dH . dH 



Qi 



dpi 



Pi 



dqt 



{'lt} = 'l,Q, {Pi}=P,P 



for both fast and slow coordinates g, Q and conjugate momenta p, P, respectively 
(hereafter, the dot over a symbol denotes time derivation). Related Hamiltonian 

Hiq,p;Q,P) = H,iQ,P)+Hj{q,p;Q) (2) 

is splitted into the slow term Hs{Q, P) and the fast one Hf{p, g; Q), latter of which 
depends on the slow coordinate also. 

In accordance with standard scheme [Hj it is naturally to pass from the fast 
variables g, p to canonical ones, being fast alternating angle (^ and slow varying 
action rf. This passage keeps invariant the first term of the Hamiltonian Q and 
transforms the second one according to relation 

H'f{^-r^,Q) = Hf{q,p-Q) + Q^-^^^ (3) 

whose explicit form is determined with generating function \l'(g; r], Q) to be defined 
by the following constrains: 

g^(9;^.Q) _ d^>{q-r^,Q) _ g^(g;r/,g) _ 

dq ^' dQ ' ^772 '^- ^^> 

Due to fast variations of the angle y?, it is naturally to consider the term (jH} 

H\r,,Q) ^ {H'f{^-r^,Q)) ^ 1- J H'f{^;7],Q)d^, (5) 



being averaged over this variations. 

To find related Hamiltonian one has to use one-valued generating function 

^^■,V,Q) = ^(g(¥^;r/,Q);r7,Q), < ^ < 27i (6) 

instead of many- valued one, \l/(g; t], Q). Then, the last factor in Eq. Q is determined 

by the relation 

(9$ d^i dq 

that expresses the chain rule. Its using gives the averaged term ^ in the following 
form: 

H'{r^,Q) = H{r^,Q) + Q/^-p^Y H{7^,Q) ^ {Hf{q,p;Q)). (8) 

As a result, usage of the canonical angle-action representation derives to the 
transformation of the averaged Hamiltonian 0: 

H,fiv;Q,P) = {H\^,r^;P,Q)) 

. «(.Q,p).Q((|i)-(4)) (.) 

where the first term 

Him Q, P) = H{r], Q) + H,{Q, P) (10) 



depends on slow varying values only. 

To find equations of motions for above slow variables one has to issue from the 
extremum condition for effective action |H] 

S,f{Q{t),P{ty,r]it)} ^ f' [P{t)Q{t) - H,f[r^{t)-Q{t),P{t))] dt (11) 

whose form is determined by the Hamiltonian Q; tj„, tf are initial and final points 
of the time. Variation of this expression with respect to the momentum derives to 
the equation of motion 

keeping initial form (PJ) (hereafter, we suppose that slow variables are vector quanti- 
ties with components Qo, Pa)- On the other hand, variation of the action (fTH) over 
the slow coordinate arrives at the equation 

that is prolonged due to an effective field with the force given by antisymmetric 

tensor 

_ dAp dAa 

and vector potential 

As usually, we use the Einstein summation rule over repeated index j3. 

Above relations pH|l - (fTH|l arrive at the conclusion of fundamental importance: 
fast variations of coordinates whose values depend on slow variables induces an 
effective gauge field [6^. For case of the Hopf bifurcation, it means that such a 
bifurcation, in the course of the system evolution, arrives at the revolving not only of 
the configuration point, but of the whole domain of the phase plane that is bounded 
by the limit cycle. In other words, the physical picture of the Hopf bifurcation means 
the phase plane behaves as a real object, but not a mathematical one. 

3 Physical picture of limit cycles 

We consider a round phase plane that is spanned on the axes of both coordinate q 
and momentum p and revolves with the angle velocity cDq and a moment of inertia /. 
From the physical point of view, the value Uq determines the frequency of external 
influence, whereas the quantity / is reduced to total action of the system under 
consideration. If the phase circle revolves as a solid plane, a phase point with 
coordinate r has the linear velocity v^ = po, r]. 

According to above consideration, the Hopf bifurcation results in the limit cycle 
creation that induces a gauge field with the vector potential fll5|) and the strength 
(HH), which are reduced to linear and angle velocities, w and u;, respectively. These 
are not equal to normal values v„ and ujq because a region bounded by the limit 



cycle revolves with different velocities due to the gauge field effect. Indeed, if one 
represents the limit cycle creation as an ordering with a complex parameter = rje^''', 
then the phase gradient v^ = sV(f, where V = d/dr, s being an elementary action, 
affects in such a manner to compensate a rotation within a domain bounded by 
the limit cycle: w = v„ — v^. In this way, the relative velocity w appears as 
a gradient prolongation: V ^ V — (i/s)w, being caused by the gauge field w. 
In opposition to the case of the solid plane revolving, ordering arrives at non-linear 
relation uj = (l/2)rot w between the angular and linear components of the revolving 
velocity. 

Well-known example of such type behaviour represents the case of revolving 
superfluid He^ [7j. Along this line, an effective potential density of the revolving 
phase plane, including a set of limit cycles, has the following form [Hj 
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E = AE{r]) + - {-isV -w)r] + -uj\ (16) 
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Within the phenomenological scheme, the density of the potential variation due to 
limit cycle creation is given by the Landau expansion 

AE{ri) = Arf + |r/^ (17) 

whose form is fixed by parameters A, B. The second term of Eq. ()16p determines 
heterogeneity energy with the gradient, being prolonged by the vector potential w 
of the gauge field. The last term is the kinetic energy of the revolving phase plane. 
Under an external influence with frequency cJq, the system behaviour is defined 
by the effective potential density 

E = E - {lujo + M) ujQ (18) 

whose value is determined with respect to the revolving plane that is characterized 
with the angular momentum M = J (cJ — uq). Steady state distributions of the order 
parameter ^^(r) and the relative velocity w(r) are given by extremum condition of 
the total value of the effective potential 

S{7]{r),w{r)} = J E(r){r), w(r))dr (19) 

where integration is fulfilled over the whole area of the phase plane. In this way, the 
boundary conditions are as follows: 

• out off a limit cycle 

r] = 0, Vr] = 0, v^ = [ujq, r], uj = ujq; (20) 

• within a limit cycle 

V = Vo, Vr] = 0, w = 0, tJ = 0; (21) 

• on a limit cycle itself 

n(-isV- w)r7 = 0. (22) 



Here, n is the unit vector being perpendicular to the hmit cycle, r/o = J —A/ B is the 
stationary value of the order parameter to be determined by the minimum condition 
of the expression (fT7|) . 

According above expressions, a disordered phase related to exterior of the limit 
cycle is characterized with the effective energy density -E(O) = —{I/2)uJq, whereas an 
ordered phase being bounded by this cycle relates to the value E^tiq) = —{\A\/2)r]Q. 
As a result, the condition -^(r^o) = E{0) of the phase equilibrium gives a character- 
istic value of the revolving velocity 




to determine an energy scale Ec = Iuj^ = \A\7]q = A^/B. Moreover, it is useful to 
introduce two lengths A, C, and their ratio n = X/^ to be determined by the following 
relations: 



A^J^. C. 



h = ^. (24) 



4|A|' ^"^2|A|' V^o' B 

Then, measuring the energy density E in units E^-, the order parameter r; - in r/o; 
the angle velocity a; - in cUc, the hnear velocities v„, w - in 2\/2XijJc-, the angular 
moment M - in 2\/2/Aa;c, and the distance r - in A, one reduces the energy density 
118j) to the simplest form 



E 



(_i^-iV - w) r/f - U - \ri') - i^o - U) ^. (25) 



Inserting this equality into the total energy ()19|) and variating the functional ob- 
tained, one finds the following equations of motion: 



k-^V't] = -[l-w^)r] + r]% (26) 

— rot rot w = r/^w. (27) 

As is known j2], the form of solutions of these equations is fixed with the pa- 
rameter K given by two last Eqs.(j21I)- ^^ usual case, the phase plane is so small 
to be realized the condition k < 2^^/^, and a single limit cycle (type of shown in 
Fig^ can be created with the form and size determined by the external frequency 
Uq. Much more rich situation is realized in the case of the so large phase plane that 
inverted condition k > 2~^/^ is fulfilled. Then, within the interval cUd < luq < uJc2, 
bounded with the limit velocities 

luK \A\ //\"\ / , , 



^c2 



= V2KuJe = \A\/s, (29) 



the mixed state is reahzed to be a set of round limit cycles periodically distributed 
over surface of the revolving phase plane. Each of these cycles has elementary 
action 27ts to reach the maximum value A^max = ^/ttC,"^ of the cycle density per unit 



of area at ujq = <jJc2- With falling down external velocity near the upper boundary 
(0 < ujc2 — cjq ^ 0Jc2) the limit cycle density decreases according to the equality 

N uoc) Vj"^ / X 

° ' (30) 



where average over the phase plane yf is connected with the revolving velocity ujq 
by the equality 

^=^^^|^('^-^o), /5 = W(^r = 0.1596. (31) 

The average value 

uJ = ujQ- lf/2n = ujo-{k- ujo)/(3{2k'^ - 1) (32) 

is smaller than external value Uq in a quantity being equal to the average of the 
plane polarization 

M=-^/2k = -{k-ujo)/P{2k^-1). (33) 

The maximum value of a revolving velocity is reached in cores of limit cycles, and 
the minimum one cUmm = ^o — V^(/^ — ^o)/(2/«^ — 1) - in the centers of triangles 
formed by cycles. The average variation of effective energy (fT^ caused by the phase 
plane revolving 



e^j^^ll^^-t\^j^l 



1 9 (k — uj) 

- + uJ - 



2 1 + i3{2k 



(34) 



is the function of the average velocity u, differentiation with respect to which results 
in Eq.dni. 

Near the lower critical value Ud, the limit cycle density A^ = {k,/2tx)uj is not so 
large and these cycles can be treated independently. Taking into account that w{r) 
varies at distances r ~ 1 and ri{r) does at r ~ k~^ <^ 1, the relative velocity are 
determined by Eq.()27|) with rf k, 1 and k ^ 1: 



w 



-K-^Kiir) (35) 



where -ft'i(r) is the Hankel function of imaginary argument. Respectively, the order 
parameter is determined by Eq. ()26|) with w = —l/nr: 

T] c^ cr at r ^ k^^, 
r/^ ~ 1 - (Kr)"^ at r > k~^ (36) 

where c is positive constant. According to Eq. ()35|1 one has w ~ —l/nr at r <^ 1 and 

w ~ —J 11 /2k? r^^/^e"*" at r ^ 1. The dependence uj{uJo) is of steadily increasing 
nature: at c^o = ojd it has the vertical tangent and with ujq growth it asymptotically 
approaches to the straight line uJ = ujq. Effective energy per one limit cycle is 
{2,11/ H?) In/t, the u) value in a cycle center is twice as large as ujci- 



4 Conclusions 

Within Hamiltonian formalism, combined consideration of both fast and slow sets 
of dynamical variables shows that averaging over the angle of the canonical pair 
angle-action induces an effective gauge field if fast coordinates depend on slow ones. 
For case of the Hopf bifurcation, it means that such a bifurcation, in the course of 
the system evolution, arrives at the revolving not only of the configuration point, 
but the whole region of the phase plane being bounded by the limit cycle. In other 
words, the physical picture of the Hopf bifurcation means the phase plane behaves 
as a real object, but not mathematical one. 

Along this line, a revolving phase plane with a set of limit cycles can be presented 
in analogy with revolving vessel containing superfiuid He^. Within framework of such 
a representation, fast varying angle is reduced to the phase y? of the complex order 
parameter cf) = rje^'^ whose module squared ■i]'^ plays a role of the action. In this 
way, a role of the vector potential of the gauge field plays the relative velocity w of 
movement of interior region of the limit cycle with respect to its exterior, whereas 
the field force is reduced to the related angle velocity u = (l/2)rot w. By this, slow 
variables are reduced to the parameters A, B of the Landau expansion ()17|1 . 

Finally, we focus upon the physical interpretation of above picture of a phase 
plane with one or more limit cycles. The usual case, when canonical pair of coordi- 
nate and momentum can be varied within the whole area of the phase plane, relates 
to the parameter k ~ 1, where the Hopf bifurcation arrives at creation of single 
limit cycle, which presents periodical variations in the self-organized system. In 
opposite case, when a domain of the coordinate and momentum variations within a 
limit cycle is much less than the phase plane area, the principle parameter is k 3> 1. 
If the self-organized system is subjected to an external periodical influence whose 
frequencies are bounded within domain Ud -h ujc2 given by Eqs. (J2H1), P^ - a set of 
non-trivial resonances is appeared. Mean values of the coordinate and momentum 
variations within these resonances relate to the centers of the limit cycles, whereas 
variations amplitude is determined by the correlation length ^ given by the second of 
definitions (j211). The first of such resonances appears at external frequency ujq = Ud 
to arrive at growth of its number to maximum value A^max ~ C~^ at the upper fre- 
quency ujq = ujc2- In this way, the resonances number varies in accordance with 
Eq.(jnni)j whereas average resonance frequency is defined by Eq. ()H2|1 which shows its 
monotonic increase with external frequency. 
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